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On  Kraichnan's  "Direct  Interaction  Approximation"  and 
Kolmogoroff's  Theory  in  Two-Dimensional  Plasma  Turbulence 
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and 

R.  N.  Sudan 

Laboratory  of  Plasma  Studies 
Cornell  University 
Ithaca,  New  York  14853 


The  nonlinear  damping  in  a  strongly  turbulent  convecting  plasma 
confuted  by  Kraichnan’s  modified  direct  interaction  approximation  and 
the  power  spectrum  are  rederived  in  a  physically  transparent  form  using 
Kolmogoroff's  theory  of  turbulence. 
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1  ? 

Recently  Sudan  and  Keskinen  ’  have  developed  a  two-dimensional 

theory  of  strong,  low- frequency'  electrostatic  turbulent  convection 

3  4 

in  a  collisional,  low  pressure  plasma  by  means  of  Kraichnan's  modified 
"direct  interaction  approximation"  (DIA) .  This  theory  explains ^  many 
of  the  features  of  the  electron  density  fluctuations  in  the  equatorial 
electrojet  observed  by  radar  scattering  techniques.^  It  is  the  purpose 
of  the  present  note  to  show  that  the  principal  results  of  the  DIA 
formalism  applied  to  the  electrojet  driven  electrostatic  turbulence 
can  be  derived  in  a  much  simpler  and  physically  intuitive  manner  from 

7 

Kolmogoroff ’ s  theory  of  turbulence. 

The  equatorial  electrojet  region  can  be  modeled  by  a  weakly  ionized 

plasma  in  a  strong  uniform  magnetic  field  Bx  subjected  to  both  a 

gradient  in  density  (3no/dz)z  and  an  electric  field  EQ  =  -  (3^/92)  z. 

The  electrons  and  ions  suffer  collisions  mainly  with  the  background 

neutrals  such  that  fl  »  u  and  D-  <  v.,  where  f2  -  and  v  .  are  the 

e  e  xx’  e,i  e,i 

respective  cyclotron  and  collisional  frequencies.  Because  of  the 
difference  in  collisionality  of  the  electrons  and  ions  the  equilibrium 

A  7 

field  Eq  gives  rise  to  a  current  density  JQy  resulting  from  EQ  x  B/B" 

7  8 

electron  drift.  Simon  and  Hoh  have  shown  that  this  configuration 
is  unstable  to  electrostatic  fluctuations  $(z)  expi(k^y-wt)  if 
Vh  .V(-e<g  >  0  which  situation  is  obtained  in  the  electrojet  during 
the  daytime. This  instability  is  analogous  to  the  Rayleigh-Taylor 
instability  by  replacing  the  gravitational  potential  with  -e<{>o.  The 
nonlinear  interaction  of  the  waves  creates  an  almost  isotropic  turbulence 
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in  the  plane  perpendicular  to  the  magnetic  field,  and  reaches  a  steady 
state  for  a  steady  electro jet  current  density  JQ  =  nee(EQ/Bo)y. 

The  electro jet  region  is  well  described  by  the  two -  fluid  equations 
applied  to  an  isothermal,  quasineutral  plasma  neglecting  both  electron 
and  ion  inertia  for  the  low  frequency  fluctuations  under  discussion. 

The  linearized  form  of  these  equations  are 


3n 


CD 


3t  ~o 


+  v  *Vn^  +  vP^»Vn  +  n  =  0, 

- n  ~e  o  o  -e  ’ 


e(E(1^  +  vp)  x  B  )  +  T TtS^/n.  *  m  v  v^  =  0, 

*  a  -  ^  o  '  o  Q  0^g 
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■ZeCECl)  +  v^  x  BQ)  +  T.Vn^V^  ♦  m^v^  =  °> 
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*  V-[nQe(Zy^  -  v^) 


V«J 


~e 


n^evj  =  0, 
^0 
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where  vQ  =  EQ/Boy  is  the  electrojet  velocity,  E^  =  -V<}»,  n^  =  Znj1^  =  n^ 
Tg,Ti  are  the  electron  ion  temperatures  and  the  rest  of  the  notation  is 
standard.  In  equilibrium  nQ(z)  =  NQ(1  +  z/L),  d4»Q/dz  =  ($q/L)  (1  +  z/L) 
and  L  is  the  electron  density  vertical  scale  height.  The  density 
fluctuation  n^  ■  Sn  expitCk^y  +  k2z)-i(wk  +  iy^jt},  with  k*B  «  0,  and 
kyL  »  l,has  a  frequency  and  damping/ growth  rate  given  by 


\  ■  k*y0/(l  + 
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where  C2  ■  (T  +  T.)/m.  and  ip  =  v  v./fl  ft. . 
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The  nonlinear  interaction  between  these  linear  modes  is  strong 

because  these  waves  are  nondispersive.  The  principal  result  of  wave 

-i  n  \ 

interaction  is  a  nonlinear  self- damping  obtained  from  DIA  theory  ’  ’  and 
expressed  by 


d2k' 


dw'  \,k-k'  Wk-k',k  Xk' ,u'fco'a),'aJk-k’+rk-k’ 


(6) 


where  the  Ws  are  matrix  coefficients  of  the  particular  model  and  the 
ensemble  average  <nf  n.  ,  ,/n2>  *  I.  5(k-k')6(<o-(o') .  This  integral 

equation  for  r.  can  be  solved,  for  known  I.  ,  by  expanding 

-  **  12 
IV  .»  about  co'-ud  »  0.  The  result  of  this  calculation  ’  is 

K”K  ,C0“U)  Kp 

given  by 
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where  k  =*  (k,0),  0  *  v./ft.  (1  +  tj>) ,  a  is  a  numerical  constant  of  order 
unity  and  the  poweT  spectrum  1^  =  f  dol^  ^  is  assumed  to  be  isotropic. 

The  basic  premise  of  Kblmogoroff ' s  mixing  length  theory  is  that 
a  particular  convective  eddy  or  mode  lasts  only  for  the  time  it  takes 
to  traverse  its  own  dimension  l,  i.e.,  the  damping  rate  is  given  by 


T  -  6v/i  ~  kfiv,  (8) 

where  6v  is  eddy  velocity.  We  express  6v  from  the  linearized  equations 
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in  terms  of  the  mode  intensity  Ij.  as  follows.  Trom  the  continuity 
equation. 


(at +  iiryo)5nk =  -iK  ■  ^0)5nk =  -iv7n0  ■  vtv (9) 


Taking  the  curl  of  the  momentum  balance  equation  (2)  we  obtain 


7  x  clvkx5o)  =  •  ?  ve  7  x  !v: 


or 


?ov’S-V  ‘  ?  ve  v  x  {,vk  -  f  Vi x  tv 


(10) 


Substituting  for  ov  from  Eq.  (5a)  and  V*6v  from  Eq.  (10)  into  Eq. 
(9)  we  obtain 
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-t— -  k*  v  uik  ~ 
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Thus,  for  kL  »  /v  we  obtain 
e  e 


fle  m  6nk  .  Vni  3nk  .  fi  6nk 
vQ  -  ve  1+*  nQ  1^  nQ  =  nQ  * 
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Now 
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and 
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rk  -  kSv  -  kvo6<CTT)2>1/2  '  Vk2li/2> 


which  is  precisely  the  result  obtained  by  previous  calculations  using 

the  modified  DIA  theory^  and  quoted  in  Eq.  (7)  except  for  a  numerical 

factor  of  order  unity.  Thus,  for  our  inodes  the  "direct  interaction 

approximation"  is  tantamount  to  accepting  the  basic  idea  behind 

7 

Kolmogorof f ’ s  theory  that  the  eddy  lifetime  is  just  the  time  it 
takes  for  the  eddy  velocity  to  cover  its  length. 

Furthermore,  the  spectrum  law  can  also  be  obtained  from  the 
elementary  Kolmogorov  argument  that  the  power  transfer  between  the 
modes  decays  only  as  the  collisional  damping,  i.e.. 


k  3E  Crk  k2V  ’  *k  kV 


2  2  2 

because  the  mode  energy  Ej.  -  ((SnjJ  -  k  1^.  Since  y^  «  -k  in  the 
damped  region  from  Eq.  (5b)  for  kL  »  1, 


!k  .  k'^  [kf 


k4/V, 


(16) 


and  kQ  is  the  cutoff  wave  number  where  Ij^  is  heavily  damped.  For 
k  «  kQ>  1^  *  k‘®/3,  since  (k/I^) dl^/dk  =  -(8/3)/[l-(k/k0)^3]  a  power 
law  is  appropriate  only  for  k  «  kQ. 

We  compare  the  results  of  our  simple  expression  with  actual  numer¬ 
ical  solutions ^  of  the  nonlinear  version  of  Eqs.  (1)  -  (4),  for  1^  in 
Fig.  1.  We  have  adjusted  our  normalization  constant  to  fit  one  point 
and  chosen  the  value  of  kQ  as  the  numerically  determined  cutoff  value. 
It  is  seen  that  the  agreement  in  shape  of  the  data  points  and  Eq.  (16) 
is  quite  good. 


In  summary,  we  have  taken  the  physical  ideas  of  KoLmogoroff  out 
of  their  context  of  three-dimensional  fluid  turbulence  and  have 
applied  them  to  a  novel  two-dimensional  situation  of  electrostatic 
plasma  turbulence  in  the  ionosphere  where  numerical  solutions 
existed.  We  have  found  that  the  ideas  carry  over  very  well  as  is 
shown  by  their  agreement  in  detail  for  the  nonlinear  self-damping  of 
the  modes  as  well  as  predicting  spectra  that  also  agree  well.  This 
indicates  that  the  mixing  length  theory  may  be  successfully  applied 
to  derive  results  for  much  more  general  situations  where  it  may  be 
difficult  or  impossible  to  carry  out  the  DIA  theory. 
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Fig. 


Figure  Captions 

1  Numerically  obtained  spectra  Ik  as  a  function  k  (Ref.  5) . 

The  continuous  line  is  obtained  from  Eq.  (16)  after  fitting 

at  one  point  marked  with  a  cross  and  the  cutoff  value  kQ  =  15. 


